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Abstract. We consider the out-of-equilibrium dynamics of the Hamiltonian Mean Field (HMF) 
model, by focusing in particular on the properties of single-particle diffusion. As we shall here 
demonstrate analytically, i/the autocorrelation of momenta in the so-called quasi-stationary states 
can be fitted by a ^'-exponential, then diffusion ought to be normal for ^ < 2, at variance with the 
interpretation of the numerical experiments proposed in Ref. ijj]. 

Keywords: Hamiltonian dynamics; anomalous diffusion; 



Several numerical studies have revealed that systems with long-range interactions 
present many interesting, and pretty unique, dynamical features which take place 
out-of-equilibrium, manifesting a strong sensitivity to initial conditions. To investigate 
these peculiar aspects and eventually develop a comprehensive theoretical framework, 
it is extremely valuable to dispose of a simple toy-model that admits an Hamiltonian 
formulation in terms of continuous variables. A paradigmatic example is represented by 
the so-called Hamiltonian Mean Field (HMF) model [2], which describes the evolution 
of A'^ particles coupled through an equally strong attractive cosine interaction. The model 
is specified by the following Hamiltonian: 



where 0,- G [0, 2n[ represents the position (angle) of the j-th particle on a unitary radius 
circle and pj stands for its conjugate momentum. The HMF model can also be seen 
as a simplification of the gravitational sheet model [HI], when considering only the first 
harmonic in the Fourier expansion of the potential. It is exactly solvable both in the 
canonical and microcanonical ensembles, leading in this case to equivalent results, and 
displays a second order phase transitions from a homogeneous to a clustered phase when 
decreasing energy or temperature. The most striking observations relate however to 
nonequilibrium dynamics. Indeed, when performing numerical simulations starting out- 
of-equilibrium, the system remains trapped in long-lived Quasi Stationary States (QSSs), 
before relaxing to thermodynamic equilibrium. In this dynamical regime the evolution is 
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INTRODUCTION 
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particularly slow and the system displays non-Gaussian momentum distributionsillH]. 
To monitor the evolution of the system, it is customary to introduce the magnetization, an 
order parameter defined as M = |M| = | ^m/|/A^, where m,- = (cos 0,-, sin0,) is the mag- 
netization vector. Simulations are often carried out for water-bag initial distributions, 
for which the one-body distribution function f{9,p,t) takes at ? = a non vanishing 
constant value only inside the rectangular phase- space domain D specified by 

D = {{d,p) E [-71, n] X [-00,00] I |0| < A0, IpI < Ap}, (2) 

where < A0 < ;r and Ap > 0. The initial magnetization Mq and the energy density U 
can be expressed as functions of AO and Ap 

_sin(A0) (Ap)2 l-(Mo)2 



This in turn implies that the initial water-bag profiles are uniquely determined by Mq 
and U, which take values in the ranges < Mo < 1 and t/ > (1 -M2) /2. 

The mean square angular displacement (7^{t) = [^KO ~ ^K^)]^ is also a quan- 
tity of interest. The scaling oc t'^ defines the diffusive behavior: 7=1 corresponds to 
normal diffusion and 7 = 2 to free particle ballistic dynamics. Intermediate cases cor- 
respond to the anomalous superdijfusive behavior. On the numerical simulations side, it 
is claimed in Ref. ^ that QSSs display anomalous diffusion with an exponent 7 in the 
range 1 .4 - 1 .5 for 0.4 < Mq < 1 

The same authors investigated the decay of the momentum autocorrelation func- 
tions and suggested that the numerical profiles are well interpolated by a so-called q- 
exponential function [5]. Based on this finding, a possible relation between the fitted 
value q and the diffusion exponent 7 has been hypothesised dll]. 

We shall here analytically prove that the aforementioned results are contradictory. 
More specifically, assuming a ^-exponential decay of the autocorrelation of momenta 
implies normal diffusion, for the range of values of q reported in Ref. [1]. 

Angular diffusion and momentum autocorrelation for the HMF model have been often 
discussed in the literature [1, @, 0, S]. In order to provide a comprehensive picture on 
this subject, we will shortly review some of these contributions in the next Section. The 
following Section is instead devoted to presenting our analytic argument. Finally we 
sum up and draw our conclusions. 



ON THE RELATION BETWEEN RELAXATION AND DIFFUSION 

FOR THE HMF MODEL 

Numerical evidence for a superdiffusive behaviour in the HMF model was first reported 
in [0]. In this paper a water-bag initial condition Q was assumed, with the particles 
positioned in di = 0, a choice which corresponds to select initially the magnetized state 



In the papers hereafter discussed, the choice U = 0.69 is always put forward, when a water-bag profile 
is initially assumed. 



Mq = 1. It was also claimed in [6] that the anomalous diffusion occurs for a transient 
out-of-equilibrium regime, when the system is trapped in a QSS. 

Successive numerical experiments [7] carried on for a homogeneous initial distribu- 
tion, i.e. Mq = 0, revealed that the anomalies in diffusion are instead associated to the 
(nonstationary) relaxation to equilibrium. The QSS should therefore be characterized by 
a normal diffusive regime: it was in fact argued that the anomalous diffusion detected in 
^ stems from finite size effects. 

Motivated by such controversial findings, the authors in [1] set down to elucidate 
the role of initial conditions in driving HMF dynamical anomalies. To this end the 
usual water-bag initial distribution was considered, with the particles confined in a 
variable portion of the unitary circle, so to span the whole interval < Mq < 1. First, the 
momentum autocorrelation function was introduced as: 

Cpit-0) = {pit)p{0))N = ^fp.WP/W, (3) 

and its decay monitored, for different values of Mq. The numerical curves were then 
interpolated with the ^-exponential: 

e,(?,T) = [l-(l-^)VT]^, (4) 

a function which arises in the realm of Tsallis' generalized thermostatistics Both T 
and q are numerically adjusted in [ll|]: q is empirically found to be 1.5 for Mq > 0.4, 
while it progressively diminuishes for smaller values of Mq. For Mq = an almost 
exponential decay is measured. In this respect, the homogeneous condition was assigned 
a special nature, and claimed to be intrinsically different from the initially magnetized 
configurations. This conclusion was further strengthened by looking at the diffusive 
properties of the system. For Mq > 0, superdiffusion is systematically observed in 
Ref. ill], while the case Mq = tends to a normal diffusive behavior when increasing the 
number of simulated particles, as previously found in [7]. Moreover, a closed analytical 
relation was conjectured between the index q and the diffusion exponent 7, based on an 
generalization of the standard Fokker- Planck equation [jofl. Such relation reads 111,190: 

7=:^ (5) 
5 — q 

and returns values of 7 similar to those measured in when the corresponding q are 
inserted. This supposed agreement is merely a coincidence, as we shall prove in the 
following. If one assumes a ^-exponential decay for the autocorrelation of momenta, 
then normal diffusion is to be exptected, for the range of q reported in [1]. Before turning 
to illustrate this point, we shall devote the remaining part of this Section to complete our 
review on the studies of diffusion for the HMF model. 

Among other interesting contributions, Ref. [10] is worth mentioning: here the con- 
troversy on the role of the initial condition is also addressed. In particular a fully mag- 
netized (Mq = 1) initial condition is investijjated. For finite sizes, a superdiffusive phase 
is indeed detected, in agreement with [ 1 , 6] . However, the local exponent 7 is shown 
to converge to the normal diffusion value 7 = 1 as the numbers of particle is increased. 



Thus, anomalous diffusion is just a finite size effect, an observation which confirms the 
scenario proposed in [7] for the homogeneous initial setting. A similar conclusion is 
also reached in [4], where the time evolution of is monitored by employing a larger 
number of particles than in previous studies. An almost normal diffusion is reported also 
for intermediate values of the initial magnetization. 

Finally, in Ref. [8] a kinetic approach that goes beyond the Vlasov approximation 
was developed. This methodology enables to esplicitly calculate the momentum auto- 
correlation function in terms of the momentum distribution at ? = 0, namely fo{p)- This 
approach applies to homogeneous initial distributions and assumes that the system is ini- 
tially in a QSS. The distribution /o(p) with algebraic tails corresponds to a correlation 
function of momenta with an algebraic decay in the long time regime, hence resulting 
in a superdiffusive behaviour. In contrast, for distributions with Gaussian tails, Cp{t;0) 
scales as \nt/t, which in turn implies an approximately normal scaling for the mean- 
square displacement of the angles. The validity of these predictions has been recently 
confirmed in IllQ. where a detailed campaign of simulations was performed. It is how- 
ever worth emphasising that the theory outlined in ^ is solely limited to homogeneous 
initial conditions. Furthermore, since the theory states that the asymptotic law of diffu- 
sion is determined by the tails of the distribution, it does not directly apply to the case 
where a water-bag initial condition for the momenta is selected. The latter has in fact no 
tails initially, although Gaussian tails rapidly develop in the course of time, as shown in 
[0], thus possibly pointing to a normal diffusion. 



^-EXPONENTIAL DECAY OF MOMENTUM 
AUTOCORRELATION IMPLIES NORMAL DIFFUSIVE 

BEHAVIOUR 

As anticipated, it is assumed in Ref. ^ that the decay of the momentum autocorrelation 
function can be fitted by a ^-exponential function These results are recurrently 
invoked to support the claim of anomalous diffusion, the exponent y and the index q 
being related through expression (|5]). As opposite to this working ansatz, we will here 
analytically demonstrate that a ^-exponential decay of the autocorrelation of momenta 
leads to normal diffusion, for the values of q that have been fitted in [1]. 

We start by noticing that the quantity o^{t) can be rewritten in terms of the correlation 
function of momenta Cp{t;w) through an exact transformation U,,^]: 

O^it) = J^dtiJ^dt2{pjih)pj{t2))N (6) 

= 2j^dsJ^ 'dwCpis-w), (7) 

where Cp{t;w) is defined as 

Cp{t-w) = {pj{t + w)pj{w))N. (8) 
Moreover, if the system is stationary, Cp{t;w) does not depend on w and hence 

Cp{t-w)=Cpit;0) (Vw>0) (9) 



then Eq.® is simplified into 

o^{t)=2 l\t-s) Cp{s;Q) ds. (10) 
Jo 

Now, let us assume a (^^-exponential decay for Cp{t;Qi), namely: 

Cp{t-Q) = [l-{\-q)t/T]-^^, (11) 
and insert the expression for Cp{s;Qi) into (flOl) 

O^it) = 2 l\t-s)\\-{l-q)s/x\^^ds (12) 
Jo 

^ ll t[l-{l~q)s/T]^'ds-2 [ s[l-{l-q)s/x]^ds. 
Jo Jo 

The first term on the right hand side of the previous equation gives: 



, t[l-{\-q)s/TY-'ids = t- / [l-z]'-^dz= (13) 

1—qJo 



= -ty^[^-t{l-q)/r]'-'+ty^ 
Z — q z — ( 



while the second yields: 



r s[l-{l-q)s/T]^^ds = -t—^[l-t{\-q)/T]^^+ (14) 
Jo 2 — q 
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■[\-t{l-q)/T]-^ + - 



{2-q){3-2q)' ' ' {2-q){3-2q) 

Finally, collecting together the two contributions, takes the form: 

The first term corresponds to a normal diffusion (7 = 1), while the latter is constant and 
can be therefore neglected for long enough times. More importantly the second term 
results in a contribution which is proportional to (1 ±Q)^ (— if ^ < 1 and + if ^ > 1), 
where v = (3 — 2q){l — q). Now, it is clear that this term causes the diffusion to be 
anomalous only if ^ > 2, otherwise resulting in a correction to the o^t behaviour. 

The functions (fT5l) have been plotted in the upper panel of Fig [T] for different values 
of q <2. In the lower panel we have reported, for the same functions, the evolution of 
the diffusion coefficient 7, which are is consistent with the results in ]^\M. 



T = no 




FIGURE 1. Upper panel: Time evolution of the angular diffusion as predicted by Eq. ( fTSl l for different 
values of q. Lower panel: Time evolution of the coefficient 7, evaluated through the logarithmic derivative 



CONCLUSION 

In this paper we have reviewed the interlaced connection between relaxational and 
diffusional properties of the HMF model. This is a debated issue which has led to highly 
controversial interpretations in the recent times. In this respect, we have here critically 
discussed the interpretation of the data presented in Ref. [1]. Our main conclusion 
goes as follows: having found a momentum autocorrelation function which is correctly 
interpolated by the relation ([TT]) . with the choice ^ ~ 1.5 [1], should necessarily lead to 
normal diffusion (asymptotic linear growth in time of the second moment). The claim 
of anomalous diffusion put forward in the same series of papers [1] seems therefore in 
contradiction with the proposed fit of the autocorrelation function. 

Indeed, such apparent contradiction could be due to finite size effects as pointed out in 
[Ivill^]. Considering that the asymptotic convergence is very slow [SIsLIio'], the estimates 
for the anomalous exponent are probably affected by the limited number of particles, an 
unfortunate fact which acts as a bias towards values of / larger than 1 . Moreover, the time 
required to converge to the normal diffusion regime is remarkably long [fllEI]. It is hence 
also possible that the time window explored by the authors in their investigations ^ is 
too short. 
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